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Abstract. We consider the error analysis of Lagrange interpolation on triangles and tetra¬ 
hedrons. For Lagrange interpolation of order one, Babuska and Aziz showed that squeezing 
a right isosceles triangle perpendicularly does not deteriorate the optimal approximation 
order. We extend their technique and result to higher-order Lagrange interpolation on both 
triangles and tetrahedrons. To this end, we make use of difference quotients of functions 
with two or three variables. Then, the error estimates on squeezed triangles and tetrahe¬ 
drons are proved by a method that is a straightforward extension of the original given by 
Babuska-Aziz. 
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I. Introduction 

Lagrange interpolation on triangles and tetrahedrons and the associated error 
estimates are important subjects in numerical analysis. In particular, they are crucial 
in the error analysis of finite element methods. Let d = 2 or 3. Throughout this 
paper, AT C M’’* denotes a triangle or tetrahedron with vertices x^, f = 1, • ■ • , d-bl. We 
always suppose that triangles and tetrahedrons are closed sets in this paper. Let \i 
be its barycentric coordinates with respect to x^. By definition, we have 0 < Ai < 1, 
— 1 - Let No be the set of nonnegative integers, and 7 = (oi, • • • , ad+i) G 
Nq'*'^ be a multi-index. Let fc be a positive integer. If jyj := J2i=i then 
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7 /A: := (ai/fc, • ■ • ,ad+i/fc) can be regarded as a barycentric coordinate in K. The 
set Yj^{K) of points on K is defined by 


{K) I |y|=fc, yeN^+l} 


Let 1 < p < oo. From Sobolev’s imbedding theorem and Morry’s inequality, we 
have the continuous imbeddings 


W^'P{K) p>d, 

W^’‘^{K) c c Vq > d, 

If d = 3, we also have the continuous imbeddings 


^,3.3/2(^) ^ ^ c Vg > 3, 

W^’P{K) C W‘^'^p/^^-p\K) C ^ C°'^-^/P{K), 1 < p < ^. 

Although Morry’s inequality may not be applied, the continuous imbedding W‘^’^{K) 
C C^{K) (d = 2,3) still holds. For the imbedding theorem, see [T], [7], and [TB]. In 
the sequel we always suppose that p is taken so that the imbedding W^^^'P{K) C 
(7° (AT) holds, that is, 


1 < p < oo, if d = 2, fc + 1 > 2 or d = 3, /c + 1 > 3, and 

3 

- < p < oo, if d = 3, fc + 1 = 2. 

Note that our discussion includes the case d = fc + l,p=l that is sometimes ignored 
in literature. 

We define the subset Tp{K) C W^^^'P{K) by 


T^{K) := {u G IF'=+i’P(iF) u(x) = 0, Vx G E'=(iF)} . 


Let Vk be the set of polynomials with two or three variables for which the degree is at 
most k. For a continuous function v G C^(K), the fcth-order Lagrange interpolation 
T^v G Vk is defined by 


From this definition, it is clear that 


Vx G s'”(AT). 


Vu G W^+^-PiK). 
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For an integer m such that 0 < m < k, is defined by 



Note that we have 



To establish the mathematical foundation of the finite element methods, we must 
show that B'^'P{K) is bounded. Many textbooks on finite element methods, such as 
those by Ciarlet [5], Brenner-Scott [B], and Ern-Guermond cni, present the following 
theorem. Let Hk be the diameter of K and pK be the radius of the inscribed ball of 


K. 


Shape-regularity. Let a > 2 be a constant. If Lk/pk < cr and hx < 1, then there 


exists a constant C = C{a) that is independent of hx such that 
Ilf; - Ixv\\i,2,K < Chx\v\2.2,K, Vn S 


The maximum of the ratio hx/px in a triangulation is called the chunkiness 
parameter [B] . The shape regularity, however, is not necessarily needed to obtain an 
error estimate for triangles and tetrahedrons. For triangles, the following estimations 
are well-known 0,0, m- 

The maximum angle condition. Let Q\, (tt/S < < tt) be a constant. If any 

angle 9 of K satisfies 9 < 9i and hx < 1, then there exists a constant C = C{9i) 
that is independent of hx such that 


Ik “ II^xV\\l, 2 ,K < Chx\v\ 2 . 2 ,K, Vu S {K). 


( 1 . 1 ) 


Later, Kffzek m introduced the semiregularity condition for triangles, which is 
equivalent to the maximum angle condition. Let Rx be the circumradius of K. 


The semiregularity condition. Letp > 1 and a > 0 be a constant. If Rx/hx < cr 
and hx < 1, then there exists a constant C = C{a) that is independent of hx such 
that 


( 1 . 2 ) 
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For tetrahedrons, the following estimation is well-known HU, [5]. 

Kfizek’s maximum angle condition. Let 02 , (tt/S < 02 < tt) be a constant. Let 
7 fc be the maximum angle of faces of a tetrahedron K and cpK be the maximum angle 
between faces of K. If jk < 02, < 02, and hx < 1, then there exists a constant 

C = C(02) that is independent of hx such that 

(1.3) \\v-l},v\\i,p,x <Chx\v\ 2 ,p,x, 'ivGW^'P{K), 2<p<oo. 

Jamet m presented a general result which covers both triangles and tetrahedrons. 
Let Ed := C be a set of unit vectors which are linearly independent. Let 

^ G be an unit vector and 0s, 0 < 0s < f be the angle between ^ and the line 
which is defined by Cs- Define 


0(Ed) '■= max min {0s}. 

jGR-ie.GBd'- 


Let iF C be d-simplex. Let N := d{d+ l)/2 and Ex be the set of N unit vectors 
that are parallel to the edges of K. Define 0x ■= niin {0{Ed)}. Note that if d = 2 

EciCEn 

and K is an obtuse triangle, then 29k is the maximum angle of K. 

Theorem 1.1 (Jamet). Let 1 < p < oo. Let m > 0, k > 1 be integers such that 
k + 1 — m > 2/p (1 < p < oo) or k — m > 1 {p = 1) if d = 2, or k + 1 — m > 3/p if 
d = 3. Then, the following estimate holds: 


(1.4) 


' ~ ^x'^\'m,p,K ^ C 


(cos 0xy 




where C depends only on k, p. 

Remark: Note that in m Theoreme 3.1] the case d = 2 and p = 1 is not mentioned 
explicitly but clearly holds for triangles. 

For further results of error estimations on “skinny elements”, readers are referred 
to the monograph by Apel [5]- 

The common idea among the above mentined estimations is that (i) show an error 
estimate for a particular type of elements, then (ii) extend it for general elements by 
affine transformation. To prove the maximum angle condition for triangles, for ex¬ 
ample, Babuska and Aziz showed the following theorem [U Lemma 2.2, Lemma 2.4]. 
Let K be the right triangle with vertices (0,0)^, (1, 0)^, and (0,1)^, and Ka be the 
right triangle with vertices (0,0)^, (1,0)^, and (0,a)^ (0 < a < 1). That is, Ka is 
obtained by squeezing K. 
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Theorem 1.2 (Babuska-Aziz). There exists a constant independent of a (0 < 
a < 1) such that < C, m = 0, 1. As an immediate consequence, we 

obtain the error estimation of Lagrange interpolation X\r on a right triangle K; for 

m = 0,1, 

|a — IxV\m,2,K < • 

Theorem O claims that squeezing a right isosceles triangle perpendicularly does 
not deteriorate the optimal approximation order of 1)^. Babuska and Aziz then 
claim that the estimate (ini) for general triangular elements is obtained by affine 
transformations. Kobayashi and Tsuchiya [T^ extended Theorem 1 1.2 1 for any p (1 < 
P < oo)- 

Now, let K denote also the reference tetrahedron with vertices (0,0, 0)^, (1, 0,0)^, 
(0,1,0)^, and (0, 0,1)^. Let Kajs be the “right” tetrahedron with vertices (0,0,0)^, 
(1,0,0)^, (0,a,0)^, and (0,0,^)'^ (0 < a,/3 < 1). 

The aim of this paper is to extend Theorem ll.2l and establish the following theorem. 

Theorem 1.3. Ifd = 2, there exists a constant Ck,m,p such that, for m = 0, • • • ,k, 

(1.5) := sup < Ck,m,p, fc > 1, 1 < p < oo. 

vCT^(Kc) \'^\k+l,p,Ka, 


If d = 3, there exists a constant Ck,m,p such that, for to = 0, • • • ,k. 


( 1 . 6 ) 

:= sup < Ck, 

vCT^iK^fj) \l’\k+l,p,Kcf3 


{ k — m = 0, 2<p< oo, 

fc = l, TO = 0, |<P< oo, 

fc>2, fc — TO>1, l<p<oo. 


Using Theorem II.31 and affine transformations, we can derive an error estimation 
on general triangles. See Section 4. 

The above mentioned estimations dTrl) . (|T 2 ]) . (fL3t . (lull) cover Theorem [L3] par¬ 
tially. We also mention that Shenk [T5] showed (II3I for p = 2, fc > 1, TO = 0, 1, and 
(11.61) for p = 2, A: > 2, TO = 0, 1. 

Because of the restrictions for to, k, and p in the above mentioned estimations, 
it seems that with fc = TO>2, l<p<2, and CH) with k = m > 2, 

1 < P < 3 have not yet been proved. To prove Theorem 11.31 we introduce the 
difference quotients of functions with two or three variables in Section [21 Then, 
Theorem 11.31 is proved in Section [3] by a method that is a straightforward extension 
of Babusla-Aziz’s original argument. The notations of functional spaces used in this 
paper are exactly the same as those in m- 
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2. Difference quotients for multi-variable functions 

In this section, we define the difference quotients for two- and three-variable func¬ 
tions. Our treatment is based on the theory of difference quotients for one-variable 
functions given in standard textbooks such as [3] and [1^. All statements in this 
section can be readily proved. 

For a positive integer k, the set C AT is defined by 

7 e 0 < |7| < kj , 

where ^jk = (ai/fc, • • • ,ad/k) is understood as the coordinate of a point in E^. 

For x..y £ E^ and a multi-index i5 £ Nq with I7I < A: — |<5|, we define the correspon¬ 
dence between nodes by 



:= x^+5 £ E^ 

For two multi-indexes 77 = (mi, • • • , rrid), 5 = (rii, • • ■ , rid), 17 < means that 

(z = 1 , • • • ,d). Also, S-rj and d! are defined by S-r] := <5! := rii! • • • n^!, 

respectively. Using A^, we define the difference quotients on E^ for / £ C^(K) by 

Let 0 := (0, ,0) £ Ng. For simplicity, we denote [xq, A^Xq] by /l^l[A^Xo]. 

The following are examples of [A'^xq]: if d = 2, 

/^[^(2’°)x(o,o)] = -^(/(x(2,0)) - 2/(X(i^o)) + /(x(0,0)))) 

0 )] = - /(x(i,o)) - /(x(o,i)) + /(x(o,o))), 

/^[A^^d)x(o g)] = ^(/(X( 2 ^ 1 )) - 2/(X(i^i)) -h /(X(o^i)) - /(X( 2 ,g)) 

+ 2/(x(i^g))-/(X(0,g))), 


and if d = 3, 

/^[A(2d.l)x((, Q 0)] = .|.(/(X(2_1_1)) - 2/(X(i,i,i)) -h /(X(g_i,i)) 
~ /(x(2,0,l)) + 2/(X(i g 1)) — /(x(g g 1)) 
~ /(X(2,l,0)) + 2/(X(i 1 g)) — /(X(g 1 g)) 
+ /(X(2,0,0)) ~ 2/(X(i g g)) -I- /(X(g g g))) 
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Let ry G Nq be such that \rj\ = 1. The difference quotients clearly satisfy the 
following recursive relations: 


If / G C^{K), the difference quotient [x^, A'^x^] is written as an integral of /. 
Setting d = 2 and S = (0, s), for example, we have 

/Mx(/.g), A(°’ 1 )x(;,,)] = k{f{xi^q+i) - /(x/,)) = d^J (^p | + dwi, 

r[x(i,p),A(°’*)xp,,)] 

/•I fWi rWs-i / I n ^ \ 

^ J J J ^ ik’k - \-Ws)] dws---dw2dwi. 

To provide a concise expression for the above integral, we introduce the s-simplex 
Sg := {(ti, ^ 2 ) ■ ■ ■ ) ^s) G \ ti > 0, 0 <ti + -- -+fs^l}) 
and the integral of g G L^{Ss) on Sg is defined by 

p pi pWi pWs — i 

/ g{wi,-■ ■ ,Wk)dWs ■■= / / ••• / g(n;i, • • • ,ri;s)dn;s • • ■dn; 2 dn;i. 

JSs Jo Jo Jo 

Let us temporarily set d = 2. Then, /®[x(;_q), A(°’®)x(/_,j)] becomes 

r [x(i.9), A( 0 ’^)xp,,)] = ^ I + ^{wi + • • • + «;,) j dWs. 

For a general multi-index (t, s), we have 
/‘+*[xp,,),A(‘>*)xp,,)] 

= J J 9)*’*)/ + —{zi + ■ ■ ■ + Zf), + • • • + Ws)^ dZtdWg. 

Let be the rectangle defined by x^ and A^Xj as the diagonal points. If i5 = {t, 0) 
or (0, s), degenerates to a segment. For v G L^{K) and with 7 = {l,q), we 
denote the integral as 

/ V ■= [ / ^ ( 7 + t(-2^i “I- Zt),Y + jiwi-\ - l-Ws) ) dZtdWg. 

Jsjst \k k k k J 
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If degenerates to a segment, the integral is understood as an integral on the 
segment. By this notation, the difference quotient is written as 


/‘+«[x^,A(‘’*)x^] 



Therefore, if u G Tp{K), then we have 


( 2 . 1 ) 



VD^*) c K. 


For the case d = 3, the integral [js v is defined in exactly the same manner. 

3. Proof of Theorem 11.31 

Let S C K he a segment. In the proof of Theorem 11.31 the continuity of the 
trace operator t defined as t : B v u|s G L^{S) is crucial. If d = 2, 

the continuity of t is standard and is mentioned in many textbooks such as [7]. For 
the case d = 3, the situation becomes a bit more complicated. If the continuous 
inclusion C C^{K) holds, the continuity of t is obvious. Even if this is 

not the case, we still have the following lemma. For the proof, see [TJ Theorem 4.12], 
[21 Lemma 2.2], and [III Theorem 2.1]. 

Lemma 3.1. Let d = 3 and S C K be an arbitrary segment. Then, the following 
trace operators are well-defined and continuous: 


t : W^'P{K) LP{S), 2<p<oo, t : W^'P{K) LP{S), l<p<oo. 


For a multi-index 6, jdj > 1, p is taken so that 


(3.1) 


2 < p < oo, if /c -I- 1 — jdj = 1, d = 3, 

1 < p < oo, if A: -I- 1 — jdj > 2, d = 3 or d = 2. 


The set C l'^l’P(iF) is then defined by 


. _ 


W'^+^~\^\’P{K) 



VD' 


ip 



By Lemma ITT] and (13.1|) . is well-defined. Note that u G Tp{K) implies d^u G 
by definition and (12.111 . 
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Lemma 3.2. We have riVk-\s\ = {0}- That is, if q G Vk-isi belongs to , 
then q = 0. 

Proof. We notice that d\mVk-\s\ = #{nfp C K]. For example, if fc = 4, d = 2, 
and |d| = 2, then dim7^2 = 6. This corresponds to the fact that, in K, there are six 
squares with size 1/4 for 5 = (1,1) and there are six horizontal segments of length 
1/2 for 5 = (2,0). All their vertices (corners and end-points) belong to E^(Ar) (see 
Figure 1). The situation is the same for d = 3. Now, suppose that v G Vk-\s\ satisfies 
/□« 9 = 0 for all □/ C K. This condition is linearly independent and determines 
<7 = 0 uniquely. □ 


Figure 1. The six squares of size 1/4 for 5 = (1,1) and the (union 
of) six segments of length 1/2 for 5 = (2,0) in K. 

The constant is defined by 


veB.l’’’ \''^\k+i-\s\,p,K 

The following lemma is an extension of [U Lemma 2.1]. 

Lemma 3.3. Let p be given by en). Then, we have < oo. 

Proof. The proof is by contradiction. Suppose that = oo. Then there exists a 
sequence C such that |wn|opx = ^ limn^oo \wn\k+i-\s\ p p = 0- 

By i Theorem 3.1.1], there exists {( 7 „} C 'Pk-\ 5 \ such that 


1 

k+l-\5\,p,K 



ll^rt + ^ 


inf llu;„ + gJI 


k+l-\S\,p,K 


qe'Pk-\5\ 


and lim \\wn + qn\\k+i-\s\ „ k = Because {iCn} C \^\’P{K) is bounded, 

n—^oo KA-± |0|,p,/v 

{qn} C Vk-\s\ is bounded as well. Hence, there exists a subsequence {qm} such that 
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converges to g G 7 ^fe-| 5 | and W^m + definition and 

the continuity of the trace operator, we have frjs Wm = 0 and 

'-^Ip 

0 = lim [ (wm +q)= [ q, Vnf C K. 

Ip Ip 

Therefore, it follows from Lemma [3.II that q = 0. This implies that 


which is a contradiction. □ 

The proof of Theorem \1.S\ The proof is an direct extension of the proof given in [H 
Lemma 2.2]. Let d = 2 initially. Define the linear transformation fa : by 


{x*,y*V = {x,ayy, G R^, 0 < a < 1, 


which squeezes the reference element K perpendicularly to Ka '■= Fa{K). Take an 
arbitrary v G and define u G by u{x,y) := v{x,ay). To 

make the formula concise, we introduce the following notation. For a multi-index 
7 = (a,b) G Nq and a real t 7 ^ 0, (a)'*'* := Let 1 < p < 00 and 1 < m < fc. 
Because u G Tp{K) and d^u G we may apply Lemma 1531 and obtain 


(3.2) 






|7|=r; 




|(5|=fe+l 


IP 

10,p,/C 


|7|=m 


0,p,K 


^ m! 


-IP 


|7|=n 


7 ! 


(fc + 1^ 

^ nUn,^VP ' ' '' 


. \r)\—k+l — m 

ml 


771 (a) 


0,p,K 


< 


|7|=m 


o,p,k 


E yW’” 

\'f\—m \ |p|=/l+ 1 —m 


^ 19 ^( 97 ^) | 7 > 


r}\ 


0,p,K 






P 

fe+1—m,p,PC 
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< 






Idyull^g 

where Ck,m,,p '■= Eaax|^|=m ^p’*- Here, we use the equality 


< CP , 


(k + iy. 

s\ 


E 


m! (fc + 1 — to)! 


7+7)=i5 

\'y\—m,\r)\—k-\-l—m 


7 ! 


77 ! 


Hence, we obtain (USD for this case. If to = 0, we have 


(3.3) 


'\o,p,K^ 


-r - 

'o.p.JC 


|5|=fc+l 


C--. Y. \3‘n\Y, 

It., 

= , ,p°'^’-^ < bI'P{KY < 00 . 


o.p.if 


E 


(fc + l)! ^gS^lP 


(5! 


0,p,K 


|u|; 

f - 

'k-\-l,p,K 


\5\^k+l 

Setting p = 00 and 1 < to < fc, we have 

iriax {(a)-^|9^w|^ 

lyl—m ^ ’ J 


(3.4) 


|'^^|m,oo,/Ca I7I— 


\''J\k+i,oo,Kc max 
|5|=fc+l 


{(«) ' l^'«lo.oo,if} 

max {(a)“''|5^w|o.oo.I?} 

|7|=m ^ ^ 


< 


niax i (a) max {(a) |5’'((9'^u)| jj} 

71=771 I 7|=/c+l —m ^ ’ ’ >» 

max {(a)“''|a'>'u|p ^ ^} 

7 =m I ’ ’ J 


< 


i?;s. {<“>■’ 1 , 1 . 1 ???-™ } 

max|,y|=^ {(«)"'>' I 

max|^l=^ {(a)-^ l^"'«lfc+i-m.oo.i?} 
max|.y|^,„ {(a)-T' 

max|,,l=^ j(a)-'y |5'’'w|o,oo.If 


< Ck', 


where Ck,m,oo '■= inax|.^|^,„ Now, the case with p = 00 and to = 0 is obvious. 
Therefore, m is proved. 
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Next, let d = 3 and repeat the above proof. Define the linear transformation 
Fa/3 : —>■ by 

{x*,y*,z*)^ = {x,ay,f3z)'^, {x,y,z)'^ 0<a,/3<l, 

which squeezes the reference tetrahedron K perpendicularly to Ka/s '■= Fap{K). 
Take an arbitrary v £ lV^^^’^(Ka/ 3 ) and define u £ by u{x,y,z) := 

v(x, ay, /3z). Let p be given by p.ll) with m = |(5|. To make formula concise, we 
introduce the following notation. For a multi-index 7 = (a, b, c) £ Ng and a real 
< yf 0, (a,/?)'’'* := Because u £ Tp{K) and d^u £ we may apply 

Lemma E31 as above. Thus, we may repeat (IX^ . (IH31) . and (EH) by replacing Ka 
with Ka/ 3 , {a)'^P with (a,/3)^^, etc. Thus, (11.61) is proved. □ 

4. Concluding Remarks 

Theorem fOl deals with only right triangles and “right” tetrahedrons. Based on 
Theorem [T31 a new error estimation of Lagrange interpolation on triangles is ob¬ 
tained in m- It should be emphasized that no geometric condition on triangles is 
imposed in Theorem 14. II 

Theorem 4.1 (Kobayashi-Tsuchiya [13]). Let K he an arbitrary triangle. Let 
1 < p < 00 , and k, m be integers such that k > 1 and 0 < m < k. Then, for the 
kth-order Lagrange interpolation on K, the following estimation holds: 

(4.1) \v - I%v\m,p,K < c h^+^-'^\v\k+yp,K = 

for any v £ W^^^''p{K), where the constant C depends only on k, p and is indepen¬ 
dent of the geometry of K. 

Any tetrahedron can be obtained from a “skinny right” tetrahedron by an affine 
transformation. To obtain an error estimate, we need to estimate the ratio of the 
maximum and minimum singular values of the Jacobian matrix of the affine trans¬ 
formation. If we obtain an expression of the ratio in terms of geometric quantities of 
the tetrahedron, a new error estimation would be obtained. The authors hope that 
they will report further development of error estimations on tetrahedrons in near 
future. 
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